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We demonstrate both experimentally and theoretically that a colloidal sphere trapped in a static 
optical tweezer does not come to equilibrium, but rather reaches a steady state in which its proba- 
bility flux traces out a toroidal vortex. This non-equilibrium behavior can be ascribed to a subtle 
bias of thermal fluctuations by non-conservative optical forces. The circulating sphere therefore 
acts as a Brownian motor. We briefly discuss ramiflcations of this effect for studies in which optical 
tweezers have been treated as potential energy wells. 



Most discussions of the dynamics of optically trapped 
particles assume at least implicitly that the forces ex- 
erted by an optical tweezer [l| are path-independent and 
therefore conserve mechanical energy. Optical forces due 
to gradients in the intensity are manifestly conservative 
in this sense 0]. Radiation pressure, by contrast, is not 
0, The experimental studies described in this Let- 
ter demonstrate that the non-conservative component of 
the optical force has measurable consequences for the dy- 
namics of optically trapped colloidal spheres. In partic- 
ular, the probability density for a sphere trapped in a 
static optical tweezer exhibits steady-state toroidal cur- 
rents, a phenomenon we call a fountain of probability. 
We use the Fokker-Planck formalism to explain how non- 
conservative forces bias random thermal fluctuations to 
induce circulating probability currents. 

I Figure [T] schematically represents the deceptively sim- 
ple system. A single colloidal sphere is drawn to the 
focus of a converging laser beam by forces arising from 
gradients in the beam's intensity [JQ. These intensity- 

\, gradient forces establish a three-dimensional potential 
energy well, V{r), determined by the local intensity, I{r). 

[ A particle trapped in this well also experiences radiation 
pressure that drives it downstream with a force propor- 
tional to /(r). In the absence of thermal fluctuations, a 
trapped particle would come to rest at a stable mechan- 

I' ical equilibrium downstream of the focus. 

Treating the displaced equilibrium point as the origin 
of an effective potential energy well is tempting but mis- 
leading. To appreciate the problem, consider a thermally 
driven trajectory such as the example shown schemati- 
cally in Fig. [TJ Were the system in thermodynamic equi- 
librium, forward and reverse trajectories around this loop 
would have equal probability. Because the light is more 
intense near the optical axis, however, radiation pres- 
sure biases the random walk in the forward direction. 
This departure from detailed balance should induce ir- 
reversible circulation in the particle's otherwise random 
fluctuations 

We demonstrate this effect by observing the motions 
of colloidal silica spheres 2.2 jim in diameter (Bangs Lab- 
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FIG. 1: (Color online) (a) Schematic representation of a col- 
loidal sphere in an optical tweezer created from a beam of light 
propagating along z. The curving arrow represents a ther- 
mally driven trajectory, (b) Measured hologram of a trapped 
sphere together with a nonlinear least-squares flt to Mie scat- 
tering theory, (c) Measured 5 min trajectory in a trap at 
Po — 7.5 mW. Emphasized portion shows 10 s. 



oratories, SS04N/7651) dispersed in water and individ- 
ually trapped in discrete optical tweezers. Each trap is 
formed from roughly 1 mW of laser light at a vacuum 
wavelength of 532 nm (Coherent Verdi) projected with 
the holographic optical trapping technique dSQ. The 
beams are brought to a diffraction-limited focus by an 
objective lens (Nikon Plan Apo 100 x, oil immersion) 
mounted in an inverted optical microscope (Nikon TE- 
2000U). We track the particle in three dimensions with 
nanometer precision using video holographic microscopy 
@i @, In-line holograms, such as the example in 

Fig. [TJb), are created by illuminating the sphere with 
a coUimated laser beam, in this case the 3 mm diameter 
beam provided by a 10 mW He-Ne laser (Uniphase) oper- 
ating at a vacuum wavelength of 633 nm. Light scattered 
by the particle interferes with the unscattered portion 
of the incident beam to produce an interference pattern 
in the microscope's focal plane. This is magnified and 
projected onto a CCD camera (NEC TI-324AII), which 
records holograms at video rates. Each frame in the video 
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FIG. 2: (Color online) (a) Effective potential localizing a col- 
loidal sphere to an optical trap at laser power Pq = 7.5 mW, 
in units of the thermal energy scale ksT. Dashed curves are 
parabolas whose curvatures yield the trap stiffness. Discrete 
points are obtained from histograms, Pj{rj), of the measured 
particle positions in 20 nm bins. The x and y results are 
offset by 1 fcsT and 0.5 fesT, respectively, for clarity, (b) 
Dependence of the effective trap stiffness on laser power. 



stream then is fit to the predictions of Lorenz-Mie scat- 
tering theory [llj to obtain each particle's position in 
three dimensions, Wt), its radius a, and its complex in- 
dex of refraction [10| . Figure [Ub) also shows the com- 
puted image resulting from fitting to the experimental 
image. A sequence of such fits yields the particle's three- 
dimensional trajectory, as plotted in Fig. [ijc). 

Were the particle to come to equilibrium within the 
trap, the probability to find it within dr of r would be 
given by the Boltzmann distribution, 



p{r)=N e^p{-pV{r)), 



(1) 



where = fcsT is the thermal energy scale at ab- 
solute temperature T and is a normalization con- 
stant. This is the basis for a popular method to cal- 
ibrate optical traps [12| , with measured particle posi- 



tions being compiled into -P(r), which then is inverted 
to yield V{r). In the further approximation that the 
trap can be modeled as a three-dimensional harmonic 
well centered at the origin, the three-dimensional proba- 
bility density can be factored into one-dimensional con- 
tributions, p{r) = Hj^i 6xp {—f3Vj{rj)), where 



1 



(2) 



The trap is then characterized by an effective spring con- 
stant kj in each of the three Cartesian coordinates. 

Figure ^a,) shows the result of applying Eq. ([2]) to 
10,000 measurements of a single particle's position in 
an optical tweezer powered by Po = 7.5 ± 1.0 mW of 
light. The difference between the trap's apparent stiffness 
= 0.467 ± 0.009 pN/^m along the axis of the hght's 
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FIG. 3: (Color online) (a) Accumulated circulation of a sin- 
gle particle's trajectory over 10,000 time steps at laser power 
Po = 7.5 mW. Dark (green) trace: experimental data. Light 
(brown) trace: simulation results for equivalent conditions. 
Inset: streamlines of the measured velocity field in the (^, z) 
plane, with the trap's center denoted by a circle, (b) Circula- 
tion rate as a function of laser power together with prediction 
of Eq. (O. 



polarization and ky — 0.400 ± 0.008 pN//im in the per- 
pendicular direction is consistent with previous reports of 
polarization effects in optical trapping 

EQEi- The 

axial stiffness, k^ = 0.080 ± 0.004 pN//im is a factor of 5 
smaller because of the comparatively weak axial intensity 
gradients [l6| . The dependence on laser power plotted in 
Fig. IHb) confirms that kj cx Pq. 

To quantify nonequilibrium effects in the particle's dy- 
namics, we consider its trajectory in cylindrical coordi- 
nates, r = (^,(/), z), centered on the trap's equilibrium 
position. During its thermally-driven exploration of the 
trap, the particle covers a roughly elliptical region in the 
(^, z) plane. A statistical bias due to radiation pressure 
along the axis should appear as a tendency of the parti- 
cle's trajectory to wind clockwise in the (^, z) plane. To 
quantify this, we define a measure of the trajectory's cir- 
culation over the r = 1/30 s the interval between video 
frames, 



_ _1_ {r{t + t) X r{t)) ■ cf, 
2vr ^(i^-iOr) {{z-{z]n 



(3) 



which is positive for clockwise circulation and negative 
for retrograde motions. The angle brackets in Eq. ([3]) 
denote an average over the trajectory. To identify trends 
in f2(i) despite large thermal fluctuations, wc further de- 
fine 



x{t) 



n{t') dt' 



(4) 



which measures the accumulation of clockwise circulation 
in the particle's trajectory. The increasing trend in x{t) 
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confirms the existence of a circulating steady state in 
which the particle completes one net cycle in roughly 9 s 
with a mean drift speed of roughly 150 nm/s. 

The measured circulation rate increases with laser 
power, as shown in Fig. Efb). The particle also becomes 
increasingly well localized as the trap stiffens, however, so 
that the range of experimentally accessible laser powers 
is limited to Pq < 15 mW by our instrumental resolution. 

Observing a convective flux in the trapped particles' 
trajectories confirms the system's departure from equi- 
librium. The nature of the nonequilibrium state is clar- 
ified by considering an idealized model of the system. 
The probability flux induced by a force F{r) acting on a 
Brownian particle is 



S{r) = fipF + DVp, 



(5) 



where p{r) is the ensemble-averaged probability density 
for finding the particle near r. We assume that the sys- 
tem reaches steady state so that p{r) does not depend on 
time. Continuity of the probability density (V • S' = 0) 
then yields the Fokker-Planck equation 



L»VV = pF ■Wp + ppW ■ F. 



(6) 



The divergence-free flux still can support circulation, 

VxS = pFx\/p + ppV X F (7) 

provided V x F ^ 0. 

We model the trap as a radially symmetric harmonic 
well with radiation pressure directed along the optical 
axis with a strength proportional to the local intensity: 



F{r) = —kr + fi cxp — 



(8) 



This model is parametrized the trap's stiffness, fc, and the 
scale of radiation pressure, /i , both of which are propor- 
tional to the laser's power, Pq. and also by the effective 
width of the trap, cr, which depends on the quality of the 
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focus and the size of the trapped particle [Ij 
For a particle whose radius a is substantially larger than 
the wavelength of light, a ^ a. Stable trapping requires 
the restoring force to exceed radiation pressure, so that 
e = fi/{ka) may be taken to be small parameter inde- 
pendent of laser power. In that case, F{r) has a point of 
stable mechanical equilibrium at z = ea. 

The lighter trace in Fig. [3lja) shows the circulation 
observed in a fourth-order Runge-Kutta Brownian dy- 
namics simulation of a particle diffusing in the force fleld 
described by Eq. ([S]), using parameters obtained from 
Fig. [21 and assuming e sa 0.1. Both the slope of x(t) 
and the magnitude of its fluctuations are consistent with 
experimental results. This agreement both supports our 
interpretation of our experimental observations and also 
indicates that Eq. ([5]) is sufficiently detailed for quanti- 
tative comparisons with our measurements. 



Expanding the probability density to first order in e, 
p{r) = pQ{r) + epi{r), yields an approximate solution [23| 
to Eq. ®, 

Po(r) = (g) ' exp (-\fikA and (9) 



pi(r) K. (3ka exp ( zpo{r). 



(10) 



Equation ^ is the result that would be obtained for a 
particle at equilibrium in a harmonic potential. Equa- 
tion (fTU)) describes the lowest-order correction due to 
spatially nonuniform radiation pressure. Because this ex- 
pression for pi{r) is obtained in the stable-trapping limit, 
/3fccr^ 3> 1, it underestimates the probability for the par- 
ticle to make large-scale excursions from the trap. The 
results that follow therefore are conservative underesti- 
mates for the infiuence of non-conservative optical forces 
on a trapped particle's dynamics. 

Even if the particle were to reach equilibrium in the 
force field F{r), the distortion of po{r) by pi{r) would 
affect measurements of colloidal forces calibrated by ther- 
mal fluctuation analysis. The first-order correction not 
only displaces the center of the probability distribution 
downstream, but also broadens it. Consequently, the ef- 
fective trap stiffness estimated from a probe particle's 
unloaded fluctuations will systematically underestimate 
the forces exerted by the trap under load, when fluctua- 
tions are suppressed. 

Our experimental results demonstrate, furthermore, 
that the distorted probability density does not come to 
equilibrium, but rather undergoes steady-state circula- 
tion. Expanding S{r) to lowest non-vanishing order in 
e, we obtain the mean circulation rate 12011 



1 ph {Pkcj^Y- 



(11) 
(12) 



This result can be generalized for an anisotropic trap 
by appropriately scaling the Cartesian coordinates in 
Eq. ®. In that case, Eq. stiU holds with the ef- 
fective spring constant set by the harmonic mean, k — 
3/(fc^^ -I- fcj^-i -I- fc~i). Given the results from Fig. H the 
measured circulation rates in Fig. [3Ib) agree well with 
the predictions of Eq. (fT2l) for e — 0.11 ± 0.04, including 

the predicted P^ dependence on laser power. 

Equation (fT^ also predicts that the circulation rate 
scales with temperature as fio ~ for (ika^ > 1. 
This confirms that the particle would not circulate at 
all were it not for thermal fluctuations. The optically 
trapped particle therefore is an exceptionally simple ex- 
ample of a Brownian motor 2l|, [22, |23[ whose ability to 



perform work relies on rectification of thermal noise. Un- 
like previous optical implementations of thermal ratchets 
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[27I [2^ , this fountain of probability requires 
no time-dependent driving, but rather is biased into mo- 
tion by the non-conservative component of the optical 
force Q. 

The fountain's efficiency as a Brownian motor can be 
estimated by computing the power dissipated into the 
water by viscous drag given the particle's mean drift ve- 
locity, v{r) = S(r)/p(r) The result 
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(13) 



depends strongly on laser power at low powers but be- 
comes independent of Pp in the strong-trapping limit. 
The maximum dissipation rate, Pmax = 5/ie^(fcBT)^/CT^, 
is a fixed proportion of the particle's thermally-driven 
velocity fluctuations. The system transduces fluctua- 
tions into heat through the bias imposed by the non- 
conservative optical force. Its efiiciency is limited by the 
material properties that determine e, and thus by the 
requirement that the particle remain trapped. 

Although this model is consistent with our measure- 
ments, other mechanisms also might cause a particle to 
circulate in an optical trap. For example, localized heat- 
ing due to absorption of light at the focal point could 
create a toroidal convection roll that would advect the 
particle. This has been demonstrated with infrared opti- 
cal tweezers operating at 1480 nm [s^l- The optical ab- 
sorption coefficient of water is five orders of magnitude 
smaller at 532 nm, however [sH. Scaling the previously 
reported [sO] convection rate accordingly suggests that 
the maximum drift speed due to thermal convection in 
our system should be no greater than 1 nm/s at the high- 
est laser powers in out study. Thermal convection there- 
fore appears unlikely to account for our observations. 

If, indeed, steady-state circulation is an intrinsic fea- 
ture of optically trapped particles' motions, the nonequi- 
librium effects we have identified should influence any 
measurement based on analysis of their thermal fluctu- 
ations. How such effects might have affected previously 
reported measurements remains to be determined. Es- 
tablishing that a colloidal particle trapped in a static op- 
tical tweezer acts as a Brownian motor also creates new 
opportunities for research in nonequilibrium statistical 
physics. For example, the nonuniform intensity in an op- 
tical tweezer also exert position-dependent torques on the 
trapped particle that we have not considered here. The 
resulting nonequilibrium behavior in the particle's rota- 
tional degrees of freedom may share features in common 
with the fountain-like translational bias we have identi- 
fied, and the two may be coupled in interesting ways. 
These considerations also may be extended to account 
for nonequilibrium effects in more general light fields. 
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